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Abstract. We estimate the effect of the experimental uncertainty in the 
measurement of the temperature of the cosmic microwave background (CMB) on the 
extraction of cosmological parameters from future CMB surveys. We find that even 
for an ideal experiment limited only by cosmic variance up to I — 2500 for both the 
temperature and polarisation measurements, the projected cosmological parameter 
errors are remarkably robust against the uncertainty of 1 mK in the firas CMB 
temperature monopole measurement. The maximum degradation in sensitivity is 
20%, for the baryon density estimate, relative to the case in which the monopole 
is known infinitely well. While this degradation is acceptable, we note that reducing 
the uncertainty in the current temperature measurement by a factor of five will bring 
it down to the per cent level. We also estimate the effect of the uncertainty in the 
dipole temperature measurement. Assuming the overall calibration of the data to 
be dominated by the dipole error of 0.2% from FIRAS, the sensitivity degradation is 
insignificant and does not exceed 10% in any parameter direction. 
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1. Introduction 

The precision measurement of the energy spectrum of cosmic microwave background 
(CMB) photons by the firas instrument [1] on board the COBE satellite [2] is one of 
the most spectacular pieces of evidence in support of the big bang theory. Owing to 
the fact that cosmological expansion does not lead to spectral distortions, but merely 
shifts the spectrum towards longer wavelengths, today's CMB spectrum still resembles 
that of a black body of temperature T [3], even though the photons have not been 
in thermal equilibrium since last scattering. As it happens, T is one of the few 
cosmological parameters that are accessible to direct measurement, without the need 
to resort to the model-dependent process of statistical inference. Combining the results 
of three independent estimation methods, the authors of reference [4] (see also [5]) find 
T = 2.725±0.001 K (at 68% c.L), an impressive accuracy AT /T of better than 0.04%. 

Since the CMB monopole T determines the present radiation density of the 
Universe, it is also a fundamental input parameter for the calculation of the temperature 
and polarisation anisotropies of the CMB. Naturally, the experimental error in T will 
also introduce a theoretical uncertainty ACfJCe in the prediction of the angular power 
spectra Ce. Depending on the scale, this uncertainty can reach a magnitude of order a 
few times 0.1% [7,8]. 

If one wants to use the observed Ce data to infer constraints on the free parameters 
of a particular cosmological model, this temperature effect ought, in principle, to be 
taken into account. In a statistically stringent Bayesian analysis, one would have to 
treat T as a free parameter and impose a suitable prior on its value instead of keeping 
it fixed. Given the statistical errors of present CMB anisotropy data [9,10], current 
parameter estimates are unlikely to be affected. However, in the near future, experiments 
such as planck [11] or cmbpol [12] will be able to measure the temperature and 
polarisation angular power spectra to an accuracy that is essentially limited by cosmic 
variance over a wide range of multipoles up to £ ~ 2500. It is therefore timely to ask 
whether parameter estimates using these high quality data sets may be compromised 
by a possibly insufficiently accurate measurement of T . 

In addition to the monopole, the analysis and interpretation of the data taken by 
full-sky CMB experiments also depends on the firas measurement of the temperature 
amplitude of the CMB dipole, r dp = 3.381 ± 0.007 mK [4]. Unlike T , the dipole 
temperature is not important for the theoretical prediction of the power spectra, but 
affects the experimental values of the s. The dipole provides a convenient way 
to calibrate detector output with the amplitude scale of fluctuations in temperature 
and polarisation. For WMAP [13] and the low- and high-frequency instruments of 
planck [14, 15], the error in the absolute calibration will be limited by the uncertainty 
in Td p , inducing a normalisation uncertainty in the angular power spectra data of 

AC cxp /c cxp = 2Ardp/Tdp „ 0. 4% [16] . 

In the present work, we determine how large an effect the uncertainty in the values 
of T and r dp will have on the estimates of cosmological parameters for the analysis of 
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future CMB data. In section 2 we outline the role played by T in the calculation of the 
anisotropies of the CMB. In section 3 we describe the technical details of our analysis, 
the results of which are presented in section 4. We summarise our results and conclude 
in section 5. A detailed account of the technicalities of generating mock CMB data is 
given in the Appendix. 



2. CMB temperature and the anisotropy spectra 

The present temperature of the cosmic microwave background is one of the basic input 
parameters for the calculation of anisotropies, affecting the evolution of the fluctuations 
during various stages of the early Universe. In particular, it determines directly the 
current photon energy density via 

P 7 ,o = ^T 4 . (2.1) 

Let us sketch briefly how the calculation of the angular power spectra will explicitly 
depend on T or p 7)0 , and to what extent these effects can mimic changes in other free 
parameters of the cosmological model. 



Baryon-to -photon ratio A fundamental input parameter in the Boltzmann equations 
for the baryon density perturbations is the baryon-to-photon density ratio, 

R^. (2.2) 

Thus, already at the level of the perturbations equations, there exists an exact 
degeneracy between T and the physical baryon density uo^ = Q^h 2 . In the tight- 
coupling limit valid before recombination, R defines the sound speed for the coupled 
baryon-photon fluid via c 2 = 1/3(1 + R), and enhances the compression phase (hence 
alternate peaks) of the acoustic oscillations [17]. The comoving sound horizon 

rv* 

r s (r?*) = / dr] c s (r)) (2.3) 
Jo 

evaluated at the time of recombination 77* governs the spacing of the acoustic peaks in 
the observed CMB anisotropies. The suppression of the anisotropy spectra at high I 
due to diffusion damping also depends explicitly on R. 

Recombination The details of the process of recombination [18,19], during which the 
photons decouple from the plasma, evidently have a significant influence on the eventual 
CMB anisotropies. Between redshifts 900 < z < 1500, the free electron fraction X c can 
be approximated by [6] 



hT l/2 
AT 

c ^ — /rP exp 



B 



(2.4) 



zT 

where h is the dimensionless Hubble parameter today, Q m is the matter density, and 
B ~ 3.9 x 10 4 K is a numerical constant. Because of the exponential dependence on the 



Effects of CMB temperature uncertainties on cosmological parameter estimation 



4 



temperature, we can expect AN e /N c ^> AT /T . From a more sophisticated calculation, 
it was shown in reference [8] that AN e /N c can be as large as 0.55% for AT /T ~ 0.04% 
within the standard ACDM model. 

Matter-radiation equality The parameter To determines not only the photon energy 
density, but also, implicitly, the neutrino energy density, and thus the total radiation 
density before the neutrinos become non-relativistic: 



Assuming the particle content of the standard model and standard neutrino decoupling, 
N e s ~ 3.046 [20]. It is apparent that any change in p r will shift the time of matter- 
radiation equality z cq , which is manifest in an enhancement especially of the first acoustic 
peak relative to the low I plateau through the early integrated Sachs- Wolfe effect. Thus, 
one can expect some degree of degeneracy between T and the physical matter density 
u> m = Q m h 2 in the CMB anisotropies. 

Projection The projection of the temperature and polarisation fluctuations onto the 
sky introduces for the observed CMB anisotropies an additional dependence on the 
angular diameter distance to the last scattering surface. For a flat geometry, 



where if (a) = lOOy/i 2 + Q m h 2 (a 3 — 1) + Q r h 2 (a 4 — 1). Evidently, projection leads to 
a degeneracy between h and T , both directly, and indirectly through u; m 's correlation 
with T . 

We conclude from this brief discussion that the parameter most degenerate with T 
is the baryon density c^b, followed by the matter density u m and the Hubble parameter 
h. Combining these three effects leads to an uncertainty in the angular power spectra 
of ACi/Ce ~ 0.2% for AT /T ~ 0.04%. For the fiducial model of section 3 we illustrate 
this uncertainty in figure 1. 

3. Methodology 

In order to answer the question whether the standard approach of keeping To fixed will 
be justified when fitting future data, we compare the results of an analysis that treats 
To and ra p as free (albeit well-constrained) parameters with those of a fit where these 
two parameters are kept constant. 

3.1. Mock data sets 




(2.5) 




(2.6) 




Following the method outlined in reference [21], we generate two sets of mock CMB 
anisotropy data, comprising the TT, TE and EE angular power spectra for multipoles 
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Figure 1. These diagrams show the difference between the angular power spectra 
under a change in To of ±lmK (thick solid lines). We also plot the effects of changing 
the baryon density Wb = ^b^ 2 by Awb/wb = 4 AT /T , holding T fixed (thin dotted 
lines). Left panel: Temperature autocorrelation. Right panel: Polarisation E-mode 
autocorrelation. 

Table 1. List of the experimental parameters assumed for the planck satellite [11]: 
$beam measures the width of the beam, At,p arc the sensitivities per pixel and v is 
the centre frequency of the channels. 



i//GHz 


^beam 


A t /(jK 


Ap/^K 


100 


9.5' 


6.8 


10.9 


143 


7.1' 


6.0 


11.4 


217 


5.0' 


13.1 


26.7 



2 < £ < 2500, assuming respectively the projected noise levels of the planck experiment 
(see table 1), and an ideal, noiseless cosmic variance limited experiment (cvl). We 
assume a sky coverage of / s k y = 0.65 in both cases. For a more detailed discussion of 
the method, we refer the reader to Appendix A. Our fiducial model is specified by the 
parameter values listed in table 2. 

In the generation as well in the subsequent analysis of the data we use the 
recombination code recfast [22,23] and ignore secondary effects such as gravitational 
lensing or the Sunyaev-Zel'Dovich effect. Let us stress that real data would require a 
less simplistic treatment of the physics of recombination (see, e.g., references [24,25]) 
and the secondary effects; failure to do so can severely bias results [26]. 

3.2. Parameter estimation 

The multi-dimensional posterior probability distributions V(9) are reconstructed using 
a modified version of CosmoMC [27], a Markov Chain Monte Carlo (MCMC) algorithm 
used in conjunction the CAMB [28] code to calculate polarisation and temperature spectra. 
For each analysis we generate eight Markov chains; their convergence is monitored using 
the Gelman and Rubin i?-parameter [29]. Our convergence criterion is R — 1 < 0.01, a 
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Table 2. In this table we show the free parameters of our model, their fiducial values 
used to generate the mock data and the prior ranges adopted in the analysis. 



Parameter Fiducial Value Prior Range 



Dark matter density 




0.104 


0.01 -> 0.99 


Baryon density 


n h h 2 


0.0223 


0.005 ^ 0.1 


Hubble parameter 


h 


0.7 


0.4^ 1 


Redshift of reionisation 




11.26 


3 -» 25 


Normalisation @ k = 0.002 Mpc -1 


ln(10 10 A s ) 


3.135 


2.7^4 


Scalar spectral index 


n s 


0.96 


0.5 -> 2 


Helium fraction 


Y Hc 


0.24 


0.1 -» 0.4 


CMB temperature 


To 


2.725 K 


see text 


CMB dipole 


Tdp 


3.381 mK 


see text 



much stricter requirement than for instance the one used by the WMAP team [30]. 

We analyse two basic models: the widely-used six-parameter "vanilla" model, and 
an extended model (vanilla+lHo) , where in addition we vary the primordial helium 
fraction. With the exception of To and Td p , we impose flat top-hat priors on the free 
parameters of the models; the limits are listed in table 2. 

For each of the models we perform the analysis with the CMB temperature and 
dipole either kept fixed at their fiducial values, or treated as free parameters. We account 
for the experimental error in To by imposing a Gaussian prior of the form 



7t(T ) oc exp 



To - 2.725 K 



(3.1) 



0.001 K 

The dipole, in principle, has no effect on the theoretical prediction, only on the data. 
However, since the absolute calibration affects polarisation and temperature data in the 
same way, on all scales, we do not need to generate new data each time Td p changes. 
Instead, we shift the calibration uncertainty to the theory side, by substituting the 
normalisation of the primordial power spectrum with 

r dp - 3.381 mK" 



1 + 2 



3.381 mK 



with a Gaussian prior on Td p , 



7r(r dp ) oc exp 



Tdp 



3.381 mKV 



(3.2) 



(3.3) 



0.007 mK 

corresponding to the result of the firas measurement. 

Note that we avoid using the popular Kosowsky-parameter 9 S [31], defined as the 
ratio of the sound horizon at recombination to the angular diameter distance to the last 
scattering surface, and fit instead the Hubble parameter h directly. Mapping between 
9 S and h as implemented in CosmoMC involves the use of a fitting formula to determine 
the recombination redshift z*, which was derived in reference [32] under a number of 
assumptions, including that of a fixed temperature, and is thus not applicable in our 
analysis. 
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Table 3. In this table we list the relative uncertainties ae/0 of the six parameters of 
the vanilla model, given in per cent. 



PLANCK PLANCKfixed 



CVL 



CVLfixcd 




0.602 
1.14 

0.859 
3.43 

0.466 

0.366 



0.607 
1.12 

0.848 
3.32 

0.465 

0.365 



0.187 

0.587 

0.387 

2.00 

0.270 

0.198 



0.149 
0.581 
0.388 
2.01 
0.274 
0.200 



ln(10 10 A s ) 



Since the future data sets considered here will be able to constrain the parameters 
of these models extremely well, one can expect the resulting posterior distribution to be 
reasonably close to a multivariate Gaussian near its mode. As a consequence, adding an 
additional parameter, such as To, with a Gaussian posterior, is unlikely to shift the point 
estimates of other parameters. Also, one would not expect the errors of uncorrelated 
parameters to be affected; only parameters that are degenerate with the new parameter 
are likely to have increased error bars. Naturally, this assumes that the Gaussian is 
actually centred around the fixed value - using a wrong value of the temperature will 
of course bias the best fit. 

The expected near-Gaussianity of the posterior also implies that different methods 
of constructing credible intervals will lead to the same results [33]. In the following, we 
will quote the standard deviation as a measure of uncertainty in parameter 9, 



where i runs over the points of the Markov chain and 9 is the mean of the 0j. This 
quantity corresponds to the width of the usual minimal 68% credible interval. 

4. Results 

4-1. Vanilla model 

The most serious potential consequence of adding extra parameters to an inference 
exercise is a shift in the parameter means, i.e., a bias in the point estimates. However, 
as expected, we find no such shift for either data set, when we compare the results from 
the fixed temperature analysis with the free temperature runs: the means differ by less 



We do find an effect on the errors though: the inferred uncertainties gq/9 of the 
vanilla model parameters are listed in Table 3. Using subsets of the full chains, we 
estimate the accuracy of these numbers to lie at the per cent level. Apart from the 
baryon density, the errors of the fixed temperature and free temperature analyses differ 




(3.4) 



than 0.1%. 
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0.0224 



g 2 0.0223 
0.0222 



2.723 2.725 2.727 

T [K] 

Figure 2. This plot illustrates the correlation between the baryon density ui^ — il^h 2 
and the CMB temperature for a fit of the vanilla model to the CVL data set. Depicted 
are the joint 2-dimensional 68%- and 95%-credible contours. 

by a few per cent at most, both for PLANCK and CVL data. This corresponds roughly to 
the expected variance of the results for multiple runs of the same model and is consistent 
with a null effect. The only significant exception is the error of the baryon density, which 
for the CVL data set is roughly 20% when taking the temperature uncertainty into 
account. The reason for this increase lies in a parameter degeneracy between To and Ut,, 
as explained in section 2 and demonstrated in figure 2. The qualitatively similar effect 
of these two parameters on the anisotropy power spectra can also be seen in figure 1. 

This mild degradation in the sensitivity to under an ideal situation indicates that 
the uncertainty in the CMB temperature measurement is, for the purpose of parameter 
estimation, sufficiently well controlled. Nonetheless, we find that a reduction in the error 
of the current temperature measurement by a factor of five will bring the degradation 
down to the per cent level. 

It is interesting to note that had we imposed instead a temperature prior of 
To = 2.726±0.005 K based on the original analysis of the firas data [34], the sensitivity 
of CVL to the baryon density would degrade by as much as a factor of 3.6 relative to 
the fixed temperature case. The degeneracies between T an d the parameters u m and 
h would also manifest as a 40% and a 10% degradation in their respective projected 
errors. Thus, through a stroke of coincidence, the current error in T of 1 mK leads 
to sensitivity degradations that are large enough to still be detectable, and yet small 
enough not to significantly limit the constraining power of even an ideal CMB survey. 
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Table 4. In this table we list the relative uncertainties ae/6 of the seven parameters 
of the extended vanilla+Yne model, given in per cent. 





CVL 


CVLfixod 


n h h 2 


0.248 


0.226 




0.574 


0.565 


h 


0.396 


0.392 




2.05 


1.99 


ln(10 10 A s ) 


0.339 


0.338 


n s 


0.317 


0.309 




1.35 


1.33 



4-2. Extended models 

While the vanilla model enjoys a large amount of popularity these days, and is generally 
used as the benchmark model for parameter estimates, it may be necessary in the 
future to consider extended models with more free parameters. One such example is 
the primordial Helium fraction, Y^e- While current CMB data are not very sensitive 
to changes in Y^ e , it will be necessary, already for PLANCK data, to include it in the 
analysis [35]. In fact, the projected sensitivity of the CMB to Yhc will rival that of 
astrophysical measurements, without being troubled by experimental systematics [36]. 

Generically, as pointed out above, adding extra parameters will tend to increase 
the uncertainties on existing parameters, provided that the data can constrain the new 
parameters well and barring unusual shapes of the posterior distribution. Our results 
for the vanilla model should thus be regarded as an estimate of the maximum possible 
effect. As can be seen from table 4, including Y^ e slightly weakens the bounds on the 
baryon density due to a degeneracy with the baryon density. As a result, the difference 
of the bounds of the fixed T and free T analyses goes down to ~ 10%. The addition 
of other parameters degenerate with Ub would further decrease the temperature effect. 

4.3. The dipole 

Equation 3.2 shows that there is a direct degeneracy between the CMB dipole and the 
inferred value of the normalisation of the initial power spectrum. From Table 3, we 
see that the relative error on the logarithm of the normalisation, ln[10 10 As], is about 
0.27% even in the most optimistic case in which the dipole is infinitely well known, 
the data is cosmic variance limited and a minimal model is assumed. This corresponds 
to a relative error in As of roughly 0.9%. Adding up this error and the dipole error 
of 0.4% quadratically, one expects an effect on the error of A s of less than 10%. This 
rough estimate is confirmed by our MCMC analysis: we find that for the minimal model 
and the CVL data set, fixing the dipole will lead one to underestimate the error of the 
normalisation A s by 8%, while the other parameters are affected by less than 1%. We 
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can thus conclude that the firas dipole measurement is sufficiently accurate for the 
purpose of future cosmological parameter inference. 

5. Conclusion 

We have shown that ignoring the uncertainty in the measurement of the present CMB 
temperature [4, 5] T = 2.725 ± 0.001 K, can affect the extraction of cosmological 
parameters from future data. However, the magnitude of this effect is rather small. 
While for projected PLANCK data it appears to be altogether negligible, one runs the 
risk of underestimating the error in the baryon density by about 20% for an ideal, 
cosmic variance limited experiment, assuming the current six-parameter vanilla model. 
For the other parameters of this model, the difference is at most at the per cent level. 
An improved measurement of To, reducing the current error by a factor of five, would 
remedy this problem. On the contrary, if one were to use the result of the original firas 
analysis [34], T = 2.726 ± 0.005 K, the effect of the temperature uncertainty would be 
much more dramatic: the projected error in the baryon density would increase by a 
factor of 3.6. Even the sensitivities to the matter density and the Hubble parameter 
would suffer some mild degradation. 

In the same vein we have also estimated the effect of the CMB dipole uncertainty. 
We found that taking into account the dipole error of 0.2% degrades the sensitivity 
to the normalisation of the primordial power spectrum by less than 10% for a cosmic 
variance limited experiment, compared to the case in which the dipole is infinitely well 
known. 

We conclude that, at least from a parameter estimation point of view, the present 
precision of CMB temperature monopole and dipole measurements is "good enough". 
However, it should be stressed that an improved measurement of the CMB spectrum 
would nevertheless be a worthwhile endeavour, for two reasons. Firstly, it offers the 
prospect for detecting possible global deviations from the blackbody spectrum, typically 
parameterised in terms of the Bose-Einstein and Compton distortions fi and y. As 
pointed out in reference [4], using state-of-the-art technology, the current 95% c.l. limits 
of < 9 x 10~ 5 and \y\ < 1.5 x 10~ 5 [3] could be improved by two orders of magnitude. 
Secondly, an actual detection of the signatures left by the process of recombination could 
serve as an additional, independent probe of cosmological parameters, such as the baryon 
density [8], as well as testing our understanding of recombination physics. 
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Appendix A. Mock data generation and the likelihood function 

We demonstrate in this section how to generate random realisations of future CMB 
data given some fiducial model, for the purpose of parameter error forecast. The 
method outlined below is essentially a generalisation of the procedure introduced in 
reference [21], and can be applied also to forecasts for, e.g., cosmic shear experiments. 

Sky maps of the CMB are usually expanded in spherical harmonics, where the 
coefficients, or the multipole moments, a^ m in the mode \i (/i = T,E,....) receive 
contributions from both the signal s^ m and the experimental noise n^ m , 

a e m — s e m + n e m - (A.l) 

Assuming the experiment has a full sky coverage and a spatially uniform Gaussian noise 
spectrum, the total covariance matrix = (a^ m a^ m ) is diagonal in the £ basis, and 
can be written as a sum of the signal Cf = {s^ m s u lm ) and noise N% u = {n^ m n v im ) power 
spectra, 

CjT = CjT + AT- (A.2) 



. Random realisations of the 

e 



Given a fiducial cosmological model 6 and the noise specifications of the experiment 
of interest, one can calculate C^\ 9q and hence Cf 
fiducial model can then be generated as follows: 

(i) Generate row vectors Ge m = {Gj m , Gj m , . . . , G™ m }, each consisting of n random 
numbers drawn from a Gaussian distribution. The number n corresponds to the 
number of observable modes (e.g., \x = T,E makes n = 2). 

(ii) The observables Ai m = {aj m , aj m , . . . , a" m } are defined as 

A£ m = Ge m L T , (A. 3) 

where L is a lower triangular matrix satisfying the relation = L ■ L T . The 

So 

components of L can be obtained from a Cholesky decomposition, so that the 



diagonal elements are given by 



1/2 



L ^=[cr a -J2Ll P \ , (A.4) 



P =i 

and the off-diagonal elements by 

1 



E VfX — — I 

Li 



0o 



L) Jj pL up , L pu = 0, (A.5) 



w \ p=l 

with v — /i + 1, /i + 2, . . . , n. 
(iii) The mock power spectra are constructed by summing the bilinear products of a^ m , 

Cr = E « m - (A.6) 

Zl -r 1 m=-i 
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To extract parameter errors from the mock data we approximate the total likelihood 
function £ as a multivariate Gaussian in the mock multipole moments a^ rn . Equivalently, 



Xeff 



-21n£ = ]T(2£ + l) 



T^C^Q) + In ■ 



n 



where we have made use of the fact that both the mock data and the noise power spectra 
are diagonal in the t basis. We approximate the effect of the mandatory sky cut near 
the galactic plane with a fudge factor / S k y , 

\Ct\ 



Xeff 



E(2^+i)/ sky 



Tr(Ci 1 C e )+ln- 



— n 



(A.8) 



where / s k y stands for the actual fraction of the sky observed after the cut. 

Finally, we note that it is also possible to perform a forecast using the fiducial 
instead of a random realisation of the fiducial model, i.e., one can set Cg equal 

in equation (A.8). This amounts to considering an average over an infinite 
number of independent realisations of the same fiducial model, and produces essentially 
similar error estimates as the more complicated procedure outlined above. 



to C, 



0o 
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